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1. Introduction 

Mean-field models, that is models in which all the particles interact with equal intensity 
regardless of distance, have been the object of considerable interest (see M for a 
review and the references therein). One of the simplest models of this kind consists 
of N particles moving on a circle coupled globally by a cosine interaction P, ||, with 
Hamiltonian: 

JV N 
i=l *J=1 

It is sometimes called the Hamiltonian Mean-Field model (HMF). Here the variables 
Pi are the momenta conjugate to di, which is the angle describing the state of the z-th 
particle. As stated above, the strength of the interparticle interaction does not depend 
on the distance and all particles interact with all others. Alternatively, one can think 
of this model as representing a mean-field approximation to the classical XY model, 
though this is certainly not a realistic model for a spin system. Its equilibrium statistical 
mechanics can be treated exactly via the Hubbard- Stratonovich transformation, as 
shown in || f|. The attractive case, corresponding to J < 0, is shown to have a 
phase transition at a certain value of the inverse temperature (3 C = 2/ J, above which 
the stable state has a clustered structure: that is, the particle density on the circle shows 
a non-trivial profile. Analogous results were also obtained using entropy maximization 
methods || |6| and, numerically, for a time discrete version of model (|1|) |7|, [|. On 
the other hand, in the repulsive case (J > 0), no phase transition is found |fl. This 
can be heuristically justified through the following physical argument. First, note that 
Hamiltonian ([]]) can be rewritten as follows: 

_/y \ / at \ 2 / ju \2 

1 v — * 2 J 



N 

I ■ 

H 




(2) 



From this follows immediately that, at least if iV is an even number, any state in 
which all spins appear in pairs (9k, tt + Ok) for all k corresponds to a ground state 
configuration of the system, all of such states having the same energy. Since the positions 
of half the angles can be chosen arbitrarily without in any way affecting the energy, 
uniform ground states will be entropically favoured, thus precluding the appearance 
of a thermodynamically stable density modulation at low temperatures, as is indeed 
rigorously established in Refs. |3|. [3J] . 

It was therefore quite intriguing to find spatially inhomogeneous solutions in the 
antiferromagnetic case [[| |^, XD | . Numerically, they have been found in molecular 



dynamics simulations starting from certain initial conditions at low temperatures. These 
so-called bicluster solutions show two paired clusters in the density, that is, particles 
appear preferentially in correlated positions (9, ir + 6). These biclusters were found to 
be extremely stable, though it is of course never easy to decide on numerical ground 
whether they are simply very long-lived or truly stationary. 
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This paper addresses the following issues: in Section |2|, we give a cursory description 
of the numerical findings concerning the bicluster. In Section |3], we show how the N 
particle system can be reduced to a system in which the particles are only linked to each 
other through the common interaction with a collective mode. We shall show that this 
phenomenon is reminiscent of the wave-particle interaction mechanism that governs a 
large amount of plasma physics phenomena O, [T^] . To make this picture even clearer 
and more quantitative, we show in Section § how this system can be further reduced to 
a one-particle time dependent Hamiltonian. Finally, in Section [|, we perform a detailed 
analysis of the one-particle model, showing that several features of the bicluster can be 
explained in these terms. Finally, in Section ^, we present our conclusions. 



2. Description of the bicluster 

In the following, we consider a system of N particles interacting according to 
Hamiltonian ([!]) with a positive value of J, corresponding to repulsion among the 
particles. We shall consider almost exclusively the following initial conditions: the initial 
velocities are all identically zero, and the angles are randomly distributed according to 
a uniform distribution on the interval [0, 2tt]. Since this is not in general an equilibrium 
position, the particles start moving and acquire typical velocities of the order of iV -1 / 2 . 
We are therefore always dealing with the low-energy dynamics of the system in the large 
N limit. Under these circumstances, it is found numerically that a bicluster forms, that 
is, the particle density becomes inhomogeneous. The early stages of this formation 
were studied analytically in Ref . |13] , using a convenient zero-temperature limit of the 



Vlasov equation. These involve the rapid formation of rather complex spatio-temporal 
inhomogeneities, so-called chevrons, which, however, disappear on a fairly rapid time 
scale. Since we are here primarily interested in the long-time behaviour, these shall not 
concern us in the following. 



Once the initial stages of growth described in Ref. [13] are over, an apparently 
stationary density p(6) arises. The time evolution of this density can be characterized 
by the following moments: 

1 - 

M k {t) = -22exp(ik6 m {t)). (3) 



m=l 

-15 



The second moment |M 2 (t)| is shown in Figure pl| for a large system (N = 10 particles). 
After some oscillations corresponding to the chevron structures described in Ref. ||13|| . 
the second moment eventually reaches a constant value of order one (though on the 
times shown in the figure, the increase is not yet quite over). Apart from quite 
rapid oscillations, which we always average over, the density p{6) settles down onto 
a provisionally stable profile, the bicluster, which is shown in Figure [B2j. 

It is found on numerical evidence |10| that the time-averaged moments in this 



quasi- stationary state are well approximated by the expression 



M (0) 
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W\ (M0) ' (4) 

whereas they are zero for odd values of k. This leads to the following analytical 
expression for the density (again after averaging over the rapid oscillations): 

Po(0) = ^-(1-M2sin0| •) (5) 
Here the moments and the density po(0) are connected by the following relation: 
r(°) _ 1 „ (a\JM. 



M"^ = —] p (6)e^d6. (6) 

The bicluster state does not last forever, at least for finite values of N. Indeed, we 
have performed additional numerical simulations of the system and could observe the 
decay of the second moment to its equilibrium value over a long time scale. Specifically, 
for iV = 50, the relaxation of |M 2 (t)| is displayed in Figure [BB] for two different energy 
densities. We see that a significant relaxation to equilibrium has taken place by a time of 
approximately 4 • 10 8 for both energies. These times rapidly increase as N increases, so 
that, in some sense, one might argue that the bicluster is some kind of equilibrium state. 
This is related to the issue of the non- commutation of the limits t —>■ oo and N — > oo in 
mean-field or Vlasov-like systems (see e.g. references [|14| , [15] for a numerical approach). 
On the other hand, the appearance of a bicluster is closely linked to the peculiar choice 
of the initial conditions: in particular, if instead of starting with zero velocities and 
random positions, one starts with random velocities of the same order as those that will 
eventually develop in the former initial condition, no bicluster formation is observed. 
Together with the fact that we are dealing with a very low-energy phenomenon in 
a system with long-range forces, this leads us to suspect that the bicluster is simply a 
phenomenon involving metastable configurations. Since it is well known that metastable 
states become infinitely long-lived in the mean-field limit, this point of view is in full 
agreement with the one involving non-commuting long-time and thermodynamic limits. 

3. A simple equivalent model 

Our purpose in this section is to map Hamiltonian (|IJ) onto a system of the following 
form: 

~ ujcpfP* NJX* 1» 

k=l 

Here P and X stand for two canonically conjugate variables. The vector (p has 
coordinates <j>k, conjugate to the pk and u(<p) is a smooth function which will be 
computed explicitly in the following. Here the <pk correspond to the angles of the original 
particles, whereas P and X correspond to an explicit separation of a collective mode. 
One therefore has the following situation: the system consists of non-interacting particles 
coupled to an oscillatory excitation, in a similar way as that familiar in the wave particle 
models in plasma physics [|ll|]. Furthermore, as we shall see, the function w((p) is an 
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average over all <pk, and therefore a slowly varying function. Under these circumstances, 
one expects a very accurate (though in general not exact) conservation of the adiabatic 
invariant related to the energy stored in the collective mode. As we shall see, if this 
invariant has an initial value much larger than its thermal average value, a bicluster 
forms and lasts as long as the invariant maintains its value. 

To show this correspondence, we start by considering system (^ as a system of 
particles in two dimensions, where all the particles are additionally constrained to move 
on the unit circle. In order to deal with the constraint, it is convenient to go to the 
Lagrangian formulation. One finds for the Lagrangian corresponding to ([!]): 

N 



t 1 j_2^ 

L = o i^( x k + yk) 



j 



2^ aK/ 2N 

k=l 

N 




+ 



+ E A ^+^-l)- (8) 

fc=l 

This has the form of a non-linearly constrained X-dimensional harmonic oscillator. We 
now separate the motion of the center of mass, since it plays a peculiar role. We therefore 
introduce: 

1 N _ 1 N _ 

k=l k=l 
X k =X k - X 1 y k = y k - X 2 . (9) 

In these new coordinates the Lagrangian reads 



1 " Mr 

l = ~ + vD + f *2 + *1 - + x l 

k=l 
N 



+ 



+ + vl + 2X i x k + 2X2y k + x\ + x 2 - 1). (io) 

*!=i 

Going back to polar coordinates and making the following approximations, which are 
valid in the limit in which the formation of the bicluster is observed, 

X 1 ,X 2 <£1 X 1 ,X 2 <1, (11) 

one finally obtains for the Lagrangian, after some manipulations shown in Appendix A, 

N X V + £ f - ^X 2 . (12) 
fe=i 

Here X and $ are related to X% and X 2 via 

X 1 = Xcos$ X 2 = Xsin$. (13) 

Let us now switch back to the Hamiltonian picture. The Hamiltonian corresponding to 
(12) is exactly given by 




JV 

2.Y7 2 2 XX 2 2 



„ u(ct>) 2 P 2 XJX 2 A 2 1^ 2 
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Here P is conjugate to X and A to $. As a final remark, we note that A is an exact 
constant of the motion, namely the total angular momentum. It can therefore always be 
set equal to zero without loss of generality. This therefore leads us to (0). The frequency 
function w(0) is given by 



V N + 2^fc=l COS 4>k 

Hamiltonian (0) is very similar in form to those developed in the theory of wave particle 
interactions in plasma physics. However, it should be emphasized at this stage that the 
system described by (0) is still an A^-dimensional system. Contrarily to the reduction 
proposed in [it] where bulk and tail (i.e. suprathermal) particles can be discriminated, 
our system is 'cold' so that all particles participate at the same level to the collective 
mode. 

A very appealing interpretation of (0) is the following. Consider a ring of mass one 
on which N beads of mass 1/N can move freely. If X is taken to be the amplitude of the 
torsional vibrations of such a ring around the y axis, and if no other motion of the ring 
is allowed, then (0) describes the motion both of the ring and the beads. Consequently, 
if we start with particles at rest and with an initial vibration entirely in the degree of 
freedom associated to the ring, we will have the ring oscillating first around its axis, 
and then driving the beads towards those regions where the ring motion is least, that 
is, towards the two points of the ring lying on the y axis. Such a structure corresponds 
exactly to what is observed for the bicluster. 

We now proceed with the calculation. Our aim is twofold. On the one hand, we 
show that equilibrium states do not contain a density modulation. On the other, we 
display a mechanism capable of leading to very long-lived inhomogeneous states. 

Let us first consider the partition function of the effective Hamiltonian (]?]). The pk 
integrals factorize and one is left with 

z = J o -al^NU^)- 1 J dPolXexp ^-^{P 2 + X 2 )^j (16) 

This, however, is easily shown not to lead to any density modulation in the angles, since 
oj(<t>) is of order one. We must therefore look beyond statistical equilibrium, something 
already indicated by the plausibility arguments mentioned in the Introduction as well 
as by exact results. 

To this end, we transform X and P locally to action angle variables, that is, we 
transform to 

1 (uj(4))P 2 NJX 2 * 



2 V NJ u(<f}) 

u(<f>)P , . 

ip = arctan NJX . (17) 

This transformation must be complemented by an approppriate transformation of the 
(f>k to make it canonical. This is computed and discussed further in Appendix B. The 
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final hamiltonian is 



I f 

k=l ^ 



uj(4>) d(f) k ) 



where the Pk are the new particle momenta after the canonical transformation. They are 
given by ( |B.3|) . Note that the interaction terms, which are linear in the pk, all average 
to zero in the fast variables ip, and the others are of order 1/N 2 . They can therefore 
be accurately taken into account via a perturbative averaging approach. The adiabatic 
theorem then states that the action I(4>) is conserved to a high degree of precision and 
over long time scales. We may therefore construct an approximately stationary state 
as follows: set an arbitrary value for the action / and, using this now as an external 
parameter, compute the resulting equilibrium state for the (pk- 

It should be noted here that recently, using a completely different approach 
involving multiple time expansion, close results have been derived by Barre et al in 
Ref. [16|. We believe that our approach complements their work by seeking to provide 



some intuition for the physical mechanisms involved in these processes. 

If we now take I to have given values in (|l"8"D and proceed to compute the Gibbs 
partition function, we can still eliminate the shifts in the P^'s and then factor the 
resulting integral out. We are then left with 

Z= I dfa ■ . . . ■ d(f) N exp [-/3w(<f>)l\ . (19) 



Jo 

Since u>((f>) is of order one, it is necessary for / to be of order N/ (3 (that is, N times larger 
than predicted by equipartition) in order to induce a density modulation of order one in 
the angles. On the other hand, if such a modulation is induced, it is straightforward to 
verify that it will indeed be a 7r-periodic modulation directed along the direction in which 
the collective mode oscillates, and therefore has most of the obvious properties observed 
in the bicluster. Strictly speaking, however, the microcanonical ensemble should be used 
in this calculation, which may possibly open a broader range of possibilities, as regions 



of negative temperature become accessible. For a discussion of these issues, see [16| . 

As stated in Section |2], various initial conditions may or may not lead to the 
formation of a bicluster. Among those that are successful, there are 

(i) Initial conditions with equispaced angles and a sinusoidal amplitude velocity 
perturbation: these initial conditions yield a bicluster always, as they put in an 
extensive (though very small) value for the action / of the collective oscillation. 

(ii) Initial conditions with random angles and zero initial velocities: For these it is 
readily verified that the actions of the collective mode are of order one, whereas the 
typical velocity of a particle decays as 1/ y/~N. From this follows that the effective 
inverse temperature (3 increases as N, thus yielding the appropriate result in (|i"9f). 
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4. Reduction to a One Particle Model 

We here show how to further reduce the model to a one-particle time dependent 
Hamiltonian. The assumptions involved are less well controlled and the results we 
obtain only coincide with the original model over fairly short times. However, in many 
features, it will turn out that the description given by this simplified model is extremely 
satisfactory. To this end, we make the following remark: in the model given, say, by 
(0), all particles are coupled to one single collective mode. The dynamics of this mode 
is determined both by an external harmonic force and by the reaction forces to all 
the particles involved. The system is therefore not a one-particle system. In order to 
reduce it to such a system, we require an approximate expression for the behaviour of 
the macroscopic mode. Once this is given, a one-particle approximation for the whole 
system immediately follows. From ([^D one obtains 

X = -oo 2 ((f))X . (20) 

Let us now consider the picture of the bicluster developed in Section |[ We had found 
there that the bicluster was a quasi-stationary state in which the action of the collective 
mode had a given (macroscopic) value and the 4> a corresponding stationary distribution. 
From this picture and ( pO"D follows that, in the bicluster, the motion of the collective 
mode can be well approximated by a simple harmonic motion. 



These observations are confirmed numerically in Figure 7 of [IIJ which shows the 
phase-points of X% in the dynamical regime where the bicluster forms. One can see that 
X\ undergoes rapid oscillations of small amplitude around the origin, with its phase 
only slowly changing in time. Note that, in reality, a slow drift in the relative phase 
between X x and X 2 is actually observed. This is not, however, due to a non-zero value 
of the angular momentum, but rather to the fact that the description of the motion as a 
two-dimensional harmonic oscillator is only approximate and not valid over large times. 

Armed with this knowledge concerning the dynamics of the collective mode, we 
can proceed to analyze the behaviour of the particles. We go back to the original 6k 
variables. The equation of motion following from (0) can now be rewritten as follows: 

9 h = J [Xi(t) sin e k - X 2 (t) cos 9 k ) . (21) 

From this, using the approximate description of the collective mode described above, 
we obtain the following one-particle approximation: 

6k = e cos uit sin 6k- (22) 

Here the small parameter e is equal to J|-X\|. Note that, by hypothesis, the amplitude 
of the collective mode X\(t) is always small, so that this last approximation is in 
fact justified (see (0)). In the next section we analyze this model in greater detail 
numerically. Before we do this, however, let us quickly show how this can be treated 
perturbatively. The dynamics ( p2]) is described by the following Hamiltonian 

H=— — e cos tut cos 6. (23) 
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Using standard averaging techniques, one finds in the limit of small e that the motion 
is well described by the following effective Hamiltonian: 

^T-s^ 0829 ' (24) 

Here is 9, from which a small rapid oscillation has been substracted and P is the 
canonicallly conjugate momentum. Explicitly 

q _ a , £ ( sin ( ut ~ °) sin(ujt + e) \ 

It follows therefore that the particle moves, at least in an average sense, in an effective 
potential of the form — cos 26, thereby leading to an attraction to the two points zero 
and 7i. 

5. Analysis of the one-particle model 

In this Section, we investigate numerically whether the dynamical features of the original 
model (0) with cold initial conditions can be accounted for by the minimal effective 
model ( p2"D with one-and-a-half degrees of freedom derived in the preceding Section. We 
first observe that we can set u — 1 without loss of generality, because this parameter 
can be absorbed in a trivial time rescaling. We then introduce the following area- 
preserving map with time-step r, that constitutes the lowest order symplectic algorithm 



that numerically solves (22) 



dn+l = O n +p n+1 T (26) 
Pn+1 =Pn~ £T COS ((n + 1) t) SU1 9 n . 



In order to reproduce the conditions of the simulations performed in [jnj, we consider 
a set of N initial conditions featuring a spatially homogeneous distribution of particles 
of zero momentum. That is, we take for all j, 1 < j < N 
2n 

^=j— andrf = 0. (27) 
We then compute the resulting iV trajectories, obtaining their superposed phase-space 



plots of figures B4 and B6. Since to = 1 the period T of the rapid oscillations of the 



collective variable X\ takes the value T = 2tt. In these simulations, we fix the time step 



t = T/100 = 27r/100 and e = 0.05. Figure |B4| shows the superposed phase-space plots 



of the N trajectories (E7|) during the initial stage from which the bicluster structure 



emerges. All trajectories wind around the elliptic fixed points 9 = 0, tc. In figure [B5 
we plot some trajectories with initial value |#o| < 7r /2, which show rapid oscillations 
of period T superposed on a much slower oscillatory motion around the minimum in 



9 = of the time averaged potential — cos 20 in Hamiltonian (|24|). The large time 
plots of figures |B6| clearly reproduce the features of the full HMF iV-degrees of freedom 
phase-space plots (e.g. figure 4 in reference |TIJ). The bicluster structure exhibits an 
overall oscillation on the period T around the zero momenta fixed points 9 = and 

9 = TV. 
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It is also possible, using the effective Hamiltonian (p4[), to get at least a qualitative 
understanding of the behaviour of the position and momentum density functions, in 
particular of their singularities. Indeed, if we assume that all particles start from rest 
and are randomly distributed on the unit circle, we find for the density p{9) 

/*2n roc 

p{9) = (27TZ)- 1 / d0 Q / alp5 [H(9,p)-H(9 ,0)] , (28) 

where Z is the microcanonical phase space volume and where we have neglected the 
difference between the transformed variables and the original ones. The integral in ( p8|) 
can be evaluated explicitly to yield 

p{9) = AT K {cos 2 9), (29) 

where K{m) is the complete elliptic integral of the first kind as a function of the 
parameter m and J\f is an appropriate normalization constant. The properties of the 
approximate solution (|29"D are very similar to the ones observed for the full iV-body 
system, e.g. we can reproduce the logarithmic singularities at 9 = 0, it. The general 



shape is also similar (see figure |B7|) , but the agreement is not as satisfactory as that of 



the empirical formula (pf), shown in figure B2 



A similar calculation for the momentum density yields 

P (p)=MK[l-^ (M<2), (30) 

where w is the appropriately scaled momentum p/(e 2 /8uj 2 ) 1 ^ 2 . Again, this has roughly 
the right form (see figure |BS| ), but, quantitatively, it is not quite satisfactory. 

Model fl2~2"|) seems therefore to be a good representation of the dynamics of the 
antiferromagnetic HMF at vanishing energy, even if the assumption of a constant 
frequency ui := 2tt/T is (slightly) violated due to self-consistency in the real system. 

Let us now show that, although simple, this model has a non trivial dynamics. 
First note that KAM theorem, which guarantees the preservation of quasiperiodic 
motions under small perturbations, cannot be applied in its classical form to this system. 
Renaming the time t as a phase variable ip, the dynamics (|22"D derives from the two- 
degrees of freedom Hamiltonian 

H (9,p, ip, u) = -p 2 + u — e cos unp cos6> (31) 

where the variable u is conjugated to time (p and does not appear in (0). KAM 
theorem requires a nondegeneracy condition on the frequency vector u) := (d p H,d u H) 
to ensure that a large set of actions (p, u) have "sufficiently irrational" frequencies. This 
condition is det (<9( P)U )Il>) ^ 0. This is trivially not satisfied here as uj = (p, 1). Yet, 
within the theory of averaging, KAM theory may be applied to an averaged Hamiltonian 



exponentially close to ( |24] ) with the result that, for e small enough, the Poincare section 
of H is filled up to a residue of exponentially small measure by invariant curves that 
are close to the level lines of H |T7], |I5| . This appears to be confirmed numerically, 
although in a rather unexpected fashion. Figure [B9| displays the enlargement of a long- 
time Poincare plot with initial conditions equally spaced on the #-axis and e = 0.05. 
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Over the time scales probed by the simulation (4 • 10 7 timesteps), the particle does not 
visit the whole phase space, but remains within a well-defined region in agreement with 
KAM predictions. In this region, however, it appears to lie on an extremely convoluted 
but smooth invariant surface. Indeed, the Poincare surface shows a dense pattern of 
parallel lines on which the points lie. These presumably represent the intersection of a 
very high order torus with the Poincare surface. Such behaviour is certainly peculiar, 
but may be well related to the degenerate nature of the unperturbed system. From this 
follows, in this particular case, that no short periodic orbits can exist in this part of the 
phase space. 

If this picture is correct, the system is not ergodic. This then strongly suggests 
that, in this reduced one-particle model, the bicluster structure lasts forever. The 
relevance of this conclusion for the N particle system, however, remains an open 
question. It is in principle conceivable that, among the class of initial conditions we 
consider, a non- vanishing measure lies on some invariant surface of the N particle 
system. However, this certainly does not follow from our previous arguments concerning 
adiabatic invariance, which only holds over very long, but not over infinite times. 
Furthermore, the simulations discussed in Section || speak against such a possibility. 

6. Conclusions 

Summarizing, we have found a physically appealing approximate description of the 
low-temperature dynamics of the repulsive Hamiltonian mean-field model in terms of 
free particles on a ring performing torsional vibrations. The physical analog of the 
ring oscillation is the collective plasmon oscillation of the variables X\ and X 2 . The 
interaction between particles and ring arises solely from the influence of the particle 
positions on the moment of inertia of the ring. This provides a straightforward 
interpretation of the clustering phenomena observed in previous work: if the initial 
condition is such that the initial ring oscillation strongly dominates the motion of the 
individual particles, then a parametric instability sets in, driving the particles towards 
the part of the ring which is at rest. This creates, as shown in Section |], an effective 
potential of the form — cos 20 in which the particles move. This same picture can also 
be used to obtain a description of the statistical mechanical equilibrium. In this case, it 
is clear that no clusters form, either in the canonical or in the microcanonical ensemble. 
Furthermore, this effective potential allows to form a qualitatively correct picture of 
both the particle density in position and in momentum space, though clear discrepancies 
remain, showing that the time-dependent nature of the problem is essential. Finally, 
we studied the chaotic properties of the time-dependent effective model. For a small, 
yet finite, modulus of the collective plasmon variable, this model was shown to present 
localization on very high order tori, and hence lack of ergodicity. This localization is 
presumably inherent to the small dimension of the system and might not occur in the 
finite iV particle model. Indeed, even if the tori present in the effective one-particle 
system actually survived in the N particle system, Arnold diffusion might still act as a 
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Appendix A. Derivation of the Approximate Lagrangian 

We start from fllQ|). We now introduce polar coordinates as follows: 



%k=Pk cos <p k y k = p k sm(j) k . (A.l) 

Here the angles <p k are measured from the vector (Xi, X2). From this follows readily via 
the cosine theorem, that 

pl + X 2 -2X Pk cos k = 1, (A.2) 

where X is defined as the norm and $ as the angle of the vector (Xi,X 2 ). Within the 
approximations stated in ([TT]) , this yields 

p k = 1 + Xcos0 fc . (A. 3) 

From this follows for the velocities 

p k = Xcos<p k - X0 fc sin0 fc = Xcos<fi k , (A.4) 

where again the last equality uses flllD and Xi,X 2 ~ 4> k . The Lagrangian fllPf ) now 
reads 

1 N N r 1 

l = ± £GS + pUD + f [* 2 + * 2 i> 2 - J* 2 } (a.5) 

k=l 

Substituting ( |A.3j ) and ([A.4Q into ( |A.5|) and using systematically the approximations 
(|TT|) yields the desired result. 

Appendix B. Transformation to Action-angle Variables 

In this appendix, we compute the canonical transformation that generates the local 
transformation to action-angle variables fllTD - We shall do this by computing its 
generating function as follows: we first determine a function Sq(I , X; (f) k ) such that 
the conditions 

f^ = P (B.l) 
dl * dX y 1 
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are equivalent to (|T7|). This is simply the generating function of the transformation 
to action-angle variables of the harmonic oscillator, where the <pk enter merely as 
parameters. If we now define 



N 



S(I, X; fa, P k ) = P k<Pk + S {I, X; 



(B.2) 



k=l 



this determines in the usual way a canonical transformation on the whole space. To 
compute it explicitly note that So depends only on the combination u(<p)I, so that 

lip dcu(4>) 



Pk 



P k + 



uj((f)) d(p k 
lip 



cos (pk sin (pk 



NJ 



COS z (pi 



(B.3) 



Note that these correction terms are of order 1/N. 
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Figure Bl. Time-dependence of |M2(t)| for a system of N — 10 5 particles at short 
times. 
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Figure B2. Bicluster density p(9) for a system of N = 10 5 particles. The full line is 
the result of numerical experiments, the dashed line is the theoretical expression (||) . 
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Figure B3. Decay of AM2(i) (|M2(i)| minus the finite N equilibrium value) for two 
systems of 50 particles with different energy densities e = H/N in log-linear scale. To 
smooth out fluctuations, AM2(t) has been averaged over time intervals of exponentially 
growing size. 
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Figure B4. Superposition of one-particle phase-space plots corresponding to the 
initial points (p7|), during the initial stage of the bicluster formation. The time unit is 
the short period T = 2n. 




Figure B6. Four snapshots at long-time within the short period T = 2ir of the phase- 
space plots of the N trajectories (|27|). The particle positions have been brought by 
periodicity to the interval [—n/2;3n/2]. 
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Figure B7. Position density p{9) at time t = 2000T of the bicluster (a) from the 
numerical integration integration of (Eq) with the 'cold' initial conditions y2fl) in plain 



line and (b) from the theoretical prediction (29) in dashed line. 




Figure B8. Density of the rescaled momentum p* — p/ y/ e 2 / (8uj 2 ) at time t = 2000T 
in the bicluster (a) from the numerical integration integration of (|2^) with the 'cold' 
initial conditions ( ^7|) in plain line and (b) from the theoretical prediction ( |30| ) in 
dashed line. 




Figure B9. Magnification of a long-time Poincare plots of system (|22| ) corresponding 
to orbits of increasing energy around the elliptic point 6 = p = 0. 



